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Raman scattering of atoms from a quasi-condensate in a perturbative regime 
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It is demonstrated that measurements of positions of atoms scattered from a quasi-condensate in 
a Raman process provide information on the temperature of the parent cloud. In particular, the 
widths of the density and second order correlation functions are sensitive to the phase fluctuations 
induced by non-zero temperature of the quasi-condensate. It is also shown how these widths evolve 
during expansion of the cloud of scattered atoms. These results are useful for planning future Raman 
scattering experiments and indicate the degree of spatial resolution of atom-position measurements 
necessary to detect the temperature dependence of the quasi-condensate. 
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I. INTRODUCTION 

Atoms scattered out of Bose-Einstein condensates can 
be an object of benchmark tests of various quantum- 
mechanical models. A prominent example is a collision of 
two counter-propagating condensates [iMzl- During the 
collision, which takes place at super-sonic velocity, atoms 
are scattered into initially empty modes, and descrip- 
tion of such process requires fully quantum treatment. 
This can be done semi-analytically in the Bogoliubov ap- 
proximation |8l-[l2] or numerically in more general cases 
|13l-[l7J . The analysis reveals strong correlations between 
the scattered atoms [3., 9-12] and sub-poissonian fluctua- 
tions of the opposite- momentum atom counts [5[ . There- 
fore, the many-body atomic states created in the colli- 
sions could have potential application for ultra-precise 
sub shot-noise atomic interferometry [18]. 

A different relevant example of atom scattering out 
of a coherent cloud takes place in a spin-1 condensate 
p^9l-[21|. In this case, a single stationary matter- wave is 
prepared in a Zeeman substate with tuf = 0. A two-body 
interaction can change the spin projection of the colliding 
pair into tuf = ±1. Recently, it has been demonstrated 
|21| that produced atomic pairs are usefully entangled 
from atom- interferometry point of view. 

Here we concentrate on another pair production pro- 
cess, namely the Raman scattering |22l, [23|. In this case, 
an ultra-cold atomic cloud is illuminated with a strong 
laser beam. As a result, an inter- atomic transition leads 
to creation of a correlated Stokes photon and atomic exci- 
tation. The scattered pairs are correlated analogously to 
those produced in the condensate- or spin-changing colli- 
sions. Raman scattering is similar to the elastic Rayleigh 
process [2J, |25| , though the Stokes photons have differ- 
ent energy then the incident light. This process has been 
widely studied theoretically [26 - 3Q| and observed exper- 
imentally in ultra-cold samples [31| and Bose-Einstein 
condensates j32l - [35| . In this work we consider a different 
source of Raman-scattered particles, namely the quasi- 
condensate, which forms in strongly elongated traps |36l - 
l38j . Due to non-zero temperature of the gas, phase fluc- 



tuations occur and they limit the spatial coherence of 
the system. This, in turn, has influence on the scatter- 
ing process. We demonstrate that one can determine the 
temperature of the parent cloud from both the density 
and the second order correlation function of the scattered 
atoms. 




FIG. 1. Scheme of atomic levels relevant in the process of 
Raman scattering. The atom, initially in state |0), absorbs 
strongly detuned pump photon with wave- vector kp. The ab- 
sorption is accompanied by spontaneous emission of a Stokes 
photon with wave- vector kg. As a result, the atom undergoes 
a transition 0^1. 



The paper is organized as follows. In Section II we 
formulate the 3-dimensional problem and introduce the 
Hamiltonian for the process of Raman scattering. We 
derive the Heisenberg equations for atoms and photons 
and introduce the relevant correlation functions. In Sec- 
tion III, basing on perturbative solution of the atomic dy- 
namics, we calculate the one-body density matrix both in 
the position and momentum representation. In Section 
IV we discuss the method for incorporating the phase 
fluctuations due to non-zero temperature of the quasi- 
condensate. Using this approach, we calculate the den- 
sity of scattered atoms both after long expansion time, 
i.e. in the far-fleld regime, and when the expansion time 
is short. Then, we turn to the second order correlation 
function. We show, how the temperature influences its 
peak height as well as the width. Some details of calcu- 



lations are presented in the Appendices. 



II. FORMULATION OF THE MODEL 

The process of Raman scattering takes places when an 
atom in a three- level Lambda configuration is illuminated 
with an intense pump beam. As a result of interaction 
with light, the atom absorbs a photon from the pump and 
undergoes an effective transition ^ 2 ^ 1 accompanied 
by spontaneous emission of a "Stokes" photon shown in 

Fig.m 

To model the phenomenon, we assume that the pump 
can be described classicahy as 



Ep{r, t) = Eo{r, t)e'(kp-r-u;pt) + ^.c, 



(1) 



where Eq is its amplitude, kp is the central wave- vector 
and Up = c|kp|. When this frequency is strongly detuned 
from the transition 0^2, the upper level can be adi- 
abatically eliminated. As a result, the process can be 
regarded as creation of a quantum of atomic excitation 
0^1 together with an emission of a Stokes photon. 

We describe the quantum state of the atoms and Stokes 
photons using two field operators 

6(k,t) = -^^e^(^— ^^)|0)«(1U, (2) 



(3) 



E^s^\r,t) = fdk e^(^-"-^l*^l^)a(k,t). 



where cjqi is the 0^1 transition frequency. The op- 
erator 6(k, t) annihilates an atomic excitation with mo- 
mentum fik, while the index a runs over all the N 
atoms in the cloud. If the majority of atoms occupy 
|0), one can apply the Holstein-Primakoff approximation 
1 , and accordingly b satisfies bosonic commutation re- 



function /i(k, k') is expressed in terms of a Fourier trans- 
form of the product of the quasi-condensate and pump 
beam fields, 

h{k,k') = hojdre-'^''+'''-''->''^{r,t)Eo{r,t), (6) 

where the coupling constant [40| is equal to 




FIG. 2. Scheme of the Raman scattering in the momentum 
space. Spontaneously emitted Stokes photon acquires mo- 
mentum ks- After many scattering events, the photons will 
form a sphere of radius ks denoted here by the dashed cir- 
cle. Due to momentum conservation, atoms scatter onto a 
sphere of radius ks as well, shifted by kp due to absorption 
of the pump photon. The width of the gray ring occupied 
by the atoms represents the uncertainty resulting from the 
momentum spread of the parent quasi-condensate. 
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(7) 



lations. Moreover, i:^+^ is the field operator of the Stokes Here, dij is the atomic dipole moment associated with 
photons. 

When a large number of atoms N occupy a single- 
particle state, one can replace summation over sepa- 
rate particles in Eq.(j2j) with an integral over the quasi- 
condensate wave- function ip{r^t). The effective Hamil- 
tonian for the process of Raman scattering is H = 
Ho + ^int, where 



^0 = [dkhwkd\k)d{k) + fdkekP{k)b{k) (4) 

is the free part, with Ck = h^k^ /2m. Also, uok = c\k\ —us 
is centered around the Stokes frequency ujs = ujp — cjqi. 
The interaction Hamiltonian 

^int = [dkdk'h{k,k') a^(k,t)6^(k^t) +H.C., (5) 

governs the desired process, where an atomic excitation 
is created together with the Stokes photon. The coupling 



the i ^ j transition and so is the dielectric constant. 
Note that in Eq. ^ we have neglected the interaction 
of the scattered atoms with the mean-field of the quasi- 
condensate. 

Below we make further, physically well justified sim- 
plifications. First, we choose the pump envelope Eo to 
be time- independent, which corresponds to a common 
situation of square pulses. Moreover, the spatial extent 
of the pump usually vastly exceeds the size of the quasi- 
condensate. Since the duration of the pump pulse is much 
shorter than the characteristic time scale of '^(r,t) dy- 
namics, the quasi-condensate wave-function can be taken 
constant and the coupling function in a frame of reference 
moving with velocity hkp/m reads 



h{k,k') = hoEotp{k^k'). 



(8) 



Here, ?/^ is a Fourier transform of wave- function ip. 

We can now derive the set of coupled Heisenberg equa- 
tions of motion for the Stokes and atomic field resulting 



from the Hamiltonian (j4]), 

ihdtd{k,t) = hwka{k,t) + fdk'h{k,k')b\k',t), (9a) 
ihdthik.t) = ekb{k,t) + fdk'h{k',k)d\k',t). (9b) 

These equations are a starting point for the analysis of 
the second order correlation function of scattered atoms, 
defined as 

G(2)(ki,k2,t) = (6t(ki,t)6t(k2,t)6(k2,t)6(ki,t)). (10) 

Since Equations © are linear and the initial state of 
scattered atoms and photons is a vacuum, then 

G(2)(ki,k2,t) = G(i)(ki,ki,t)GW(k2,k2,t) 

+ |GW(ki,k2,t)p, (11) 

is a function of the (one-body) density matrix, which 
reads 



G(i)(ki,k2,t) = (5t(ki,t)S(k2,i)). 



(12) 



Its diagonal part n(k, t) = ^^-^^(k, k, t) represents the 
momentum distribution of scattered atoms. 

In the following section we derive analytical expres- 
sions for the density matrix in momentum and position 
representations treating the atom-photon interactions in 
the perturbative manner. 



III. PERTURBATIVE SOLUTIONS FOR ATOMS 

When the number of scattered atom-photon pairs is 
small, one can solve Eq. (po]) perturbatively in the cou- 
pling constant ho defined in Eq. ([7]), 

ihdtU^\k,t) = ekb^^\k,t) -^fdk'h{k',k)d\k')e'^^'\ 

(13) 
where a^k') = a^k'jO). As we argue in Appendix El 
since Ck ^ ficok, the first order solution can be written as 



U^'>{k,t)=e-'^b{k) 



(14) 



te * 2ft 
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.^ dk'h{k, k')d^ (kO sine ( ^ j e'^ . 

This expression is used to calculate the first order corre- 
lation function ([12]) of scattered atoms. 

The measurement of positions of scattered atoms is 
performed as follows. First, the initial wave-packet of 
the quasi-condensate interacts with the pump beam for a 
time tp and atoms scatter out of the mother cloud. Then, 
the system freely expands for a time t/, and finally the 
positions of atoms are recorded. If tf is sufficiently long 
and the system reaches the far-field regime, positions of 

atoms r^ are related to their wave- vectors k^ by k^ = xr^ • 

* ^ -^ ^ htf 

In the present work we compare the two possible ex- 
perimental situations, when the system either is or is not 




FIG. 3. Graphical representation of the density integral ([16]) 
for k ^ ksGz. The gray ring represents the halo of scattered 
atoms, in analogy to Fig. [J] When k = /cse^, as in the main 
part of the picture, the integration samples the majority of 
the atomic cloud, leading to value of the peak height, (a) 
When |k| < ks, the integration samples the tails of the quasi- 
condensate and thus the density of scattered atoms does not 
vanish, (b) However, when |k| > kg, the tails do not con- 
tribute to the integral anymore and the density drops rapidly 
with growing k. 



in the far field. In the former case, as we argued above, 
it is sufficient to calculate the density matrix ([T2]) in the 
momentum space just after the interaction ends. In the 
latter, we provide an expression for G*^^^ as a function of 
expansion time in position space. 



A. Momentum-dependent density matrix 

In order to calculate the momentum-dependent den- 
sity matrix, note that for typical interaction times, the 
"sine" function appearing in Eq. (fH|) . which is peaked 
around k' = ks^ has much smaller width than the Fourier 
transform of the condensate function, which, via Eq. ([8]), 
enters h{k^k'). Therefore, one can fix the length of the 
photon wave- vector to be equal to kg. Using the defi- 
nition from Eq. (fT2]) and the solution from Eq. ([T4|) we 
obtain the density matrix in the momentum representa- 
tion, 

G^i^ki, k2) = ajdO.' ^*(ki + /c,nO^(k2 + /c^n^, (15) 



where we omitted an irrelevant phase factor and a = 
-^^|/ioPl^oP^s- The integration is performed over all 
the directions of the unit vector n^ Since the above per- 
turbative expression, apart from a trivial scaling of a 
with tp, is time-independent, we have skipped the time 
argument of G^^\ All the intermediate steps leading to 
the above solution are presented in Appendix [Bl 

By setting ki = k2 = k we obtain the density of the 
scattered atoms 



/' 



n(k) =a dfl \^{k + ksn)\' 



(16) 



which is directly related to the momentum distribution 
of the quasi-condensate. Integration over all directions 
of Stokes photon momentum nks is characteristic for a 
spontaneous regime, where photons scatter isotropically. 
Note also that since the quasi-condensate wave-function 
from Eq. ([8]) is expressed in the reference frame moving 
with the velocity h'kp/m, in the laboratory frame scat- 
tered atoms form a halo centered around k^ vector , see 
Fig.H 
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FIG. 4. Density of scattered atoms for k ~ ks^z for differ- 
ent temperatures and normalized by the value of the peak 
height. Each curve is an average over 400 realizations of the 
phase noise. The black solid line is calculated for the quasi- 
condensate at T = O.lnK, the dotted blue line for T =200nK 
and the dashed red line for T =960nK. 



IV. NUMERICAL RESULTS 

In this section we calculate the correlation functions 
([15]) and (p^S]) using realistic experimental parameters. 
First, we briefly describe a numerical method for simu- 
lating phase fluctuations present in a strongly elongated 
ultra-cold bosonic gas. 
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FIG. 5. Density of scattered atoms for ^ ~ /csGz for different 
temperatures and normalized by the value of the peak height. 
The time of flight is t/ = 300 ms. Each curve is an average 
over 400 realizations of the phase noise. The black solid line is 
calculated for the quasi-condensate at T = O.lnK, the dotted 
blue line for dashed T =200nK and the dashed red line for 
T =960nK. 



B. Position-dependent density matrix 

To deal with situations when tj is not sufficiently long 
for the system to enter the far-field regime we provide 
an expression for the density matrix in position space 
as a function of the expansion time t/, which up to an 
irrelevant phase factor reads 

GW(ri,r2,t/) 



(27r)6 



(17) 
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Here, the functional S is given by 



d{Sk) e-'^^'f^Pik + Sk). (18) 



^(k) 



The derivation of Eq. (p!8|) is presented in detail in the 
Appendix O Note that ([18]) resembles ([15]) except that 
ijj is replaced with Slip]. The atomic density is obtained 
by setting ri = r2 = r in Eq. ([T8]) and reads 

2 



i{r,tf) 



(2^)6 
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htf 



k.n 



(19) 



We also underline that for sufficiently long t/, 6'['^(k)] 
'?/i(k), so the far field is reached. 



A. Quasi-condensate 

We apply the above model to the process of Raman 
scattering of atoms from A^ = 10^ metastable ^He bosons 
with the atomic mass m = 6.65 x 10~^^kg and the scat- 
tering length a = 7.5 x 10~^ m. The atoms are confined 
in a harmonic potential with the radial and the axial fre- 
quencies equal to cj^ = 27r x 1500^, cj^ = 27r x 7.5^. 
Such an elongated gas is called a "quasi-condensate" due 
to presence of the phase fluctuations along the z axis 
[36 --38]. 

To account for the quasi-condensate fluctuations, we 
use the method introduced in [36|, 133]. First, we evaluate 
numerically the density profile of the pure condensate by 
finding a ground state of the stationary Gross-Pitaevskii 
equation 



ficp{r) 



where g 



2m 



mujz 



mujz 



z^^g\^{r)\' 



(20) 



47rh^ 



- and /i is the chemical potential. 



Next, we construct the quasi-condensate wave function 
?/^(r) by imprinting a phase (j){z) onto the pure condensate 



function, '0(r) = |(/?(r)|e*'^^^\ where 




g^^(i + 2)(2i + 3) ^(1.1) 



'-tf 



«i- (21) 



Here, 



hujz 



VJiJ + 3) is ^h^ energy of the low-lying 



.(1,1) 



axial excitations and P ' is the Jacobi polynomial. 
Also, z^f is the axial density width given by Thomas- 
Fermi approximation. 
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FIG. 6. The peak of the second order correlation function 
^^^'*(0) as a function of temperature. Different curves cor- 
respond to averaging over various number of realizations. 
Clearly, g^^'(^^ > 2 and the inequality is saturated only for 



(jj. :^ ClO/cg. All the numerical results presented below 
are obtained by calculating the relevant physical quan- 
tity for a single realization of the phase noise (j){z) and 
then averaging over many such realizations. 

We can now estimate the free expansion time t/, at 
which the system enters the far-field regime in the z 
direction. The velocity spread of the quasi-condensate 
/S.kz is approximately 2 mm/s, while the initial size is 
2ztf — 1 mm. Therefore, the far-field condition would be 
tf > 2zti//^kz = 0.5 s. 

Let us now comment on the consistency of the above 
approximations. Equation (|2T]) was derived in [37| under 
assumption, that the quasi-condensate has a Thomas- 
Fermi density profile in all three dimensions, valid when 
the non-linear term dominates the Gross-Pitaevskii equa- 
tion. To simplify our calculations, we model radial wave- 
functions with Gaussians, which is true in the quasi- ID 
limit, when the radial trapping potential dominates over 
the non-linear term. However, when we plot the numer- 
ically evaluated ground state of Eq. (|2Q]) . it turns out 
to be in an intermediate regime, and could be equally 
well modeled with either a Gaussian or a Thomas-Fermi 
shape. Therefore, the above method can be regarded 
as an appropriate approach in such a transitional case. 
Since our main goal is to demonstrate the general behav- 
ior of the density and the correlation functions of scat- 
tered atoms as a function of T, we believe that this ap- 
proximate method is sufficiently precise for the purpose. 



B. Density of scattered atoms as a function of T 



The phase fluctuations result from randomness of aj , 
which is drawn from a Gaussian distribution with zero 
mean and the variance given by occupation of the j-th 
mode 



(«?) 



1 



exp{ej/kBT) 



(22) 



Consequently, the temperature of the gas enters the dy- 
namics of the Raman scattering. Due to presence of the 
^-dependent phase factor ([2T]) in the quasi-condensate 
wave-function, the momentum distribution along the z 
axis broadens as the temperature grows. 

The quasi-condensate is illuminated with an intense 
laser beam with the wave- vector equal to kp = 5.80 x 
10^ m~^. Since the Stokes photon and the pump wave- 
vectors are similar, we set ks ^ kp. 

Our final simplification regards the form of the con- 
densate density profile. A simple numerical check shows 
that the ground state of Eq. (|2Q|) can be approximated 
by a Gaussian function in the radial direction, so that 



^(r) 



Me-4(^'+y').^(,)e^ 



(23) 



The axial function (p{z) is found numerically by set- 
ting X = y = in Eq. (|2Q|) and the Gaussian fit gives 



First, we present the numerical results for the momen- 
tum distribution of scattered atoms, as given by Eq. ([16]). 
We use the reference frame co-moving with a velocity 
fikp/m, hence according to Fig. [2] the density is centered 
around k = with the radius equal to kg . We investigate 
the momentum density as a function of kz in a vicinity 
of k = /cgG^. This quantity, via Eq. (fT6|) . samples the z- 
dependence of the momentum distribution of the quasi- 
condensate and therefore may provide some information 
on its temperature. 

In Fig. [3] we schematically show which ksU vectors con- 
tribute to the integral Eq. (p!6|) with k c^ kg^z- When 
kz = kg, as in the main part of the Figure, the integra- 
tion runs approximately through the center of the cloud. 
When kz < ks^ as in the inset (a), the tails of the quasi- 
condensate still contribute to the density so the integral 
does not vanish rapidly as we move away from kz = kg. 
On the other hand, when kz > kg-, shown in (b), the in- 
tergal does not sample the tails anymore and the density 
of scattered atoms quickly drops with growing kz . 

This simple graphical interpretation is readily con- 
firmed in Fig. m where we present the result of numerical 
integration of Eq. (fT6|) with the quasi-condensate func- 
tion obtained for various temperatures and averaged over 
400 realizations of the phase noise. As expected, for 
T = 0.1 nK the density is peaked around kz — kg, and 



is largely extended for kz < /Cg, while it drops immedi- 
ately as kz > kg. As the temperature grows, the den- 
sity widens substantially due to increased width of the 
quasi-condensate in the momentum space. Therefore, the 
density of scattered atoms, when measured along the kz 
axis, could be used to determine the temperature of the 
mother cloud. 

Next, we investigate the dependence of the atomic den- 
sity measured in position space after a typical expansion 
time of tf = 300 ms [3|-[5|. At this time, as argued in the 
previous section, the system has not yet entered the far- 
field regime. We set ^ o:^ ks^z in Eq. (p!9]) to make a di- 
rect comparison with previous results and evaluate the in- 
tegral numerically. We observe that when the expansion 
time is finite, contrary to the far-field regime considered 
above, the density is less sensitive to the temperature. 
Also, we notice that for high T results in both cases are 
more similar. This is because at large temperatures, the 
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FIG. 7. The dependence of the second order correlation func- 
tion g^'^\Ak) on the length of the wave- vector A/c in a vicin- 
ity of ksez. The black solid line is calculated for T = 0.1 nK, 
the dotted blue for T = 200 nK while the dashed red for 
T = 960 nK. The number of realizations was 400 for each 
curve. 



momentum distribution of the quasi-condensate broad- 
ens and so the far-field condition is satisfied at earlier 
times. 



C. Dependence of correlation functions on the 
temperature 

We now investigate the impact of the phase fiuctua- 
tions on the correlations of scattered atoms. We begin 
with the far-field expression ([15]) and using Eq. (pT]) we 
calculate the normalized second order correlation func- 
tion 



<7^'Hki,k2) = 



(n(ki)),(n(k2))/ 



(24) 



where (•)0 denotes averaging over many realizations of 
the phase <j){z). In order to be consistent with the re- 
sults of the previous section, we concentrate on a re- 
gion of wave- vectors close to /c^e^. Namely, we set ki = 

ks^z + ^^z and k2 = kgBz — ^ez and analyze g^'^\/S.k). 
First, note that for A/c = 0, according to Eq. (pT]) we have 
^(2)/Q) ^ 2(n^^ ^^^^^ ^ ^ G^^\ksez,ksez). Since the 

variance (n^)^ — (n)^ is non-negative, then ^*^^^(0) > 2. 
The inequality is saturated only in the absence of noise 
fiuctuations. To picture the impact of the temperature 
on the height of the peak of the second order correlation 
function, in Fig. [6] we plot ^'^^^(0) for various tempera- 
tures. We clearly notice the change of the height of the 
peak as soon as T > 0. 
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FIG. 8. The dependence of the second order correlation 
function g^'^\/\z^tf) in the position space in a vicinity of 
^- — ksBz as a function of Az. The black solid line is cal- 
culated for T = 0.1 nK, the dotted blue for T = 200 nK while 
the dashed red for T — 960 nK. The number of realizations 
was 400 for each curve. 

Next, in Fig.[7]we plot g^'^\Ak) as a function of Ak for 
various temperatures averaged for 400 realizations. We 
observe that apart from the change of the peak height 
^•^^^(0), the wings of the correlation function broaden, due 
to increased momentum width of the quasi-condensate at 
higher temperatures. When T c^ 0, the correlation func- 
tion oscillates in the momentum space. This behavior is 
determined by a Fourier transform of the Thomas- Fermi 
profile in the kz direction. At larger temperatures, the 
phase fiuctuations smear out the fringes. 

We now switch to the finite expansion time regime and 
evaluate the normalized second order correlation function 
in position space. Namely, we use the definition (pT]) , 
where the G^^^ function is calculated using Eq. (p!8|) . giv- 
ing 



^(2)(ri,r2,t/) 






(25) 



In analogy to the far-field case, we set ri 



- fvaGz 



-Gz and r2 = 



m 



— kg^z 



-e^. In Fig. [8] we plot 



g^'^\Az^tf) as a function of Az for three different val- financial support of the National Science Centre, 
ues of temperature T and tf = 300 ms. Equivalent ly to 
the results obtained in the momentum space, we have 
g^'^\0^tf) > 2 and agian we observe that this equal- 
ity is saturated only for T c^ 0. However, differently 
from the previous case, the oscillations of the correlation 
function at low temperatures are not present, since the 
Thomas- Fermi profile is a smooth function of z. Also, 
the broadeding of the correlation function is much less 
pronounced. 



Appendix A: Perturbative solution for the atoms 

In this Appendix we present details of the derivation 
of Eq. (p!4|) . First, we introduce a solution of Eq. (|9a|) in 
absence of coupling, i.e. 



a(k, t) = a(k)e' 



-iuJkt 



(Al) 



V. CONCLUSIONS 

We have analyzed in detail the properties of the field 
of atoms scattered out of a quasi-condensate in the Ra- 
man process. We have demonstrated that the den- 
sity of scattered atoms, when measured in the far-field 
regime, strongly depends on the temperature of the 
quasi-condensate. However, this dependence is much 
weaker, when the expansion time is finite. Furthermore, 
we have calculated the second order correlation function 
in both expansion time regimes. In each case, ^^^^ broad- 
ens with growing T, although in the latter the effect is less 
pronounced. The presence of the temperature-induced 
phase fluctuations can be also deduced from the peak 
height of ^^^^. While for the pure condensate, ^^^^(0) = 2, 
this value can substantially rise at higher T. 

In summary, the measurements of the position of scat- 
tered atoms could provide some information on the tem- 
perature of the mother quasi-condensate. Nevertheless, 
physical quantities such as the density or the correlation 
functions do not change drastically in a wide range of 
temperatures T G [0,1] fiK. If the experiment is aimed at 
determining the temperature of the quasi-condensate, it 
requires very high spatial resolution low detection noise 
and long expansion time. 

Note that in our calculations we have neglected the 
possible impact of the atomic transition rules on the field 
of scattered atoms. In the case of particular atomic tran- 
sitions, due to polarization of the Stokes field, some scat- 
tering directions are forbidden. We underline, that such 
effect could be easily taken into account by modifying the 
coupling function /i(ki,k2). 
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This expression, when inserted into Eq. ()9b|) . gives a first 
order equation of motion for the atomic field, which reads 

5,/3(k,t) = i fdk'h{k\k)d^{k')e'^^^^^^^\ (A2) 

where b{k,t) = p{k,t)e~'^^^^^^. Integration over time 
gives 

/3(k,t) = ^(k,0) + ^ldkh{k,\^)a\k!) X 

xsinc(^(K, + |)0e^(-^' + *)i (A3) 

The typical values of the kinetic energy of scattered 
atoms are Ck — ^7^, while the photon energies are of 
the order oiujk' — hkpC. Since mc/h ~ 1.9- lO^^m"^ is of 
the order of the inverse of "Compton wavelength" of an 
atom one can drop the dependence on the atom energy in 
the "sine" function. We now express the above equation 
in terms of operator b and arrive at Eq. (p!4|) . 



Appendix B: C*^^^ in momentum space 

In this Appendix we derive Eq. (p!5]) . Using the pertur- 
bative solution from Eq. (p!4|) and the definition of G^^^ 
from Eq. (p!2|) . we obtain up to an irrelevant phase factor 



Ukdp 



G^^\ki,k2) = -^ /dk'/i*(ki,kO/i(k2,kOsinc' 

(Bl) 
Typically the duration of the pump pulse is of the order of 
5/is. For such value, ^~ 1.5-10"^— is much smaller 

f^ ' ctp m 

than the width of the coupling function /i, which, via 
Eq. ([8|) , is related to the Fourier transform of the quasi- 
condensate function. Since Uk' is centered around cjg, we 
can set |k'| = ks in the coupling function and perform 
the integral over /c^ This way, we obtain 

G(i)(ki,k2) ^^^Idn h*{kuksn')h{k2,ksn'). 

(B2) 
where k' = kgU.'^ \n'\ = 1 and JdQ' denotes integration 
over a solid angle pointed by n^ Using the definition of 
the coupling function from Eq. (|8|) we arrive at Eq. (p!5]) . 



Appendix C: G^^^ in position space 



where a 



27rtpfc . T^ |2|7^|2 



ch'^ 



Eo\'\ho\' and 



In this Appendix, we calculate the first order correla- 
tion function in the position space after tf time of the 
free expansion. We employ a reasonable approximation, 
that the pump duration time tp is much shorter than 
tf. In this case, the position-dependent correlation func- 
tion G(ki, k2, tj) is simply given by the following Fourier 
transform of Eq. (|B2)) 



^(r,n,t/) 



(27r)3 



e '-^ 



Changing the variables to k 



G^ ^(ri,r2,tj) = a I dVt ^*(ri, n, tj)^(r2, n, tj), (CI) up to an irrelevant phase factor. 



V;(k + fc,n). (C2) 



(5k gives Eq. ([T8|), 
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